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Abstract 

Background: Tlie purpose of tliis study was to evaluate tine role of ventilatory parameters in maximal exercise 
performance in elite soccer players. 

Methods: From September 2009 to December 2012, 90 elite soccer players underwent evaluation of lung function 
test and ergospirometry by means of an incremental symptom-limited treadmill test. Results were analyzed according 
to i) maximal exercise velocity performed (Hi-M: high-performers, > 18.65 km/h; Lo-M: low-performers, < 18.65 km/h) 
and ii) usual role in the team. 

Results: Hi-M showed higher peak minute ventilation {VEpeak- 158.3± 19.5 vs 148.0± 18.54 L/min, p = 0.0203), and 
forced expiratory volume at first second (5.28 ± 0.50 vs 4.89 ± 0.52 liters, p < 0.001) than Lo-M, independently of 
playing role. Moreover, a significant correlation between peak oxygen uptake and VE (r = 0.57, p < 0.001) was found. 

Conclusions: Ventilatory response plays a role in the assessment of exercise capacity in elite soccer players. 
Keywords: Exercise testing. Field performance. Team sports. Training, Ventilatory parameters 



Background 

Playing soccer is the result of combination of technical, 
tactical, psychological and athletic skills. To assess soccer 
player performance model, studies have used match ana- 
lysis systems [1-5], and laboratory and field tests [6-9]. 
Other studies addressed different roles [10,11], different 
data tests [12,13], and relationships between endurance 
field tests and match analysis data [14]. 

Soccer is a sport characterized by more than 1,000 
unpredictable and acyclic changes in activity, each oc- 
curring every 3 to 5 seconds, involving up to 40 sprints, 
tackles and jumps per match [15]. Decelerations, kicking, 
dribbling, and tackling are additional actions required as 
well [16]. Therefore, the physical effort imposed on the 
players is elicited by all these gestures: as a consequence 
soccer is a highly physiologically demanding sport, with 
additional stress resulting by frequent matches and high 
load training sessions are quite high [17-21]. 
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Computerized systems for time-motion analysis have 
shown that during a match elite soccer players are able 
to perform 2 to 3 km of high-intensity running (at a speed 
>15 km/h) and about 0.6 km sprinting (>20 km/h): these 
distances are 28% and 58% respectively greater than those 
covered by intermediate level professional players [10]. 
Technical and tactical skills in soccer are highly dependent 
on the player physical capacity [22,23]. 

Relevant effort has been devoted to the evaluation of 
muscular strength, power, maximal speed, aerobic and 
anaerobic endurance, but poor attention has been ad- 
dressed to ventilatory capacity. In most studies cardio- 
respiratory response to maximal exercise was evaluated 
using peak oxygen uptake {V02peak) and the correspond- 
ing anaerobic threshold value [23-31]. To the best of 
our knowledge, currently available studies do not focus 
on ventilatory function as evaluation tool of professional 
soccer players performance. 

In the present study we hypothesized that ventilatory 
capacity can be a determinant of soccer players exercise 
capacity. In order to demonstrate our hypothesis the 
relationship between ventilatory parameters and exer- 
cise tolerance was evaluated in a group of elite Italian 
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soccer players using an incremental symptom-limited 
cardiopulmonary exercise test on a treadmill Several 
different measurements reflecting cardiovascular, respira- 
tory, and metabolic response can be recorded during the 
exercise test. The assessment of ventilatory parameters 
can contribute to improve training programmes and, con- 
sequently, athletic performance. 

Methods 

Subjects 

From 2009 to 2012, ninety professional soccer players from 
five Italian serie A soccer teams were evaluated in the 
period September-December of each year after completing 
pre-seasonal training programme in the frame of routine 
medical evaluations. Each subject provided informed signed 
consent to the use of their data for scientific purposes. 

As expected, no player reported either smoking habit or 
any relevant disease, with negative chest physical examin- 
ation. Players were unaware of the aim of the study and 
researchers performing analysis of results were blind to 
players' identity. 

Procedures 
Lung function test 

Dynamic lung volumes were assessed by means of a 
pneumotachograph (V-Max Encore, Yorba Linda, CA, 
USA). Predicted values were those of American Thoracic 
Society [32] . Maximal Voluntary Ventilation (MVV) was 
estimated multiplying Forced Expiratory Volume at first 
second (FEVi) value by 40 [33]. 

Electrocardiography 

Resting and exercise electrocardiography (EKG) was as- 
sessed in upright position by means of a 10-lead electro- 
cardiograph (Cardiosoft, GE medical systems, Fairfield, CT, 
USA) applied on the cardiac screening (six precordial leads) 
and the posterior wall of the chest (four peripheral leads). 

Exercise test 

An incremental symptom-limited exercise test was 
performed on a treadmill (Runrace 900, Technogym, 
Gambettola, Italy) under EKG and pulse oximetry moni- 
toring. Subjects standing on the treadmill breathed 
through a mask. A continuous "ramp" protocol at con- 
stant grade (1%) (starting from 8 km/h, increasing 
speed by 1 km/h every 60 seconds) was used. The test 
was stopped when subjects complained of exhaustion. 
Exercise tolerance was evaluated as the maximal speed 
reached (Maximal Exercise Velocity: MEV), adjusted ac- 
cording to a modified Kuipers equation (Equation 1) [34]. 



MEV = V/ + 




where v/ represents the speed achieved at the last exer- 
cise step and n the number of seconds attained during 
the last stage. 

Gas measurements 

The following variables were measured at peak exercise 
through breath by breath analysis of inhaled and exhaled 
gases, by mass flow meter and fast- responding gas ana- 
lyzer (V-Max Encore, Yorba Linda, CA, USA): Oxygen 
uptake {VO2) and VO2 normalized to body weight, and 
its relationship with heart rate (HR) (pulse oxygen or 
VO2/HR), CO2 production (FCO2), the physiological 
dead space to the tidal volume ratio (V^/V^), minute 
ventilation {VE), maximal VE {VEyyiax) expressed as 
the highest VE value recorded either during exercise or 
at the first time recovery phase, breathing respiratory 
reserve (BRR% expressed as VE^ax to MVV ratio). The 
anaerobic threshold (AT) was estimated by V slope 
method and ventilatory equivalent method [35] . Anaerobic 
phase time (APT) was defined as the time spent during 
VC02lV02> 1. Predicted values of HR were computed 
according to Tanaka et al. [36] . 

Statistical analysis 

Subjects were categorized into four groups according to 
their role in the team as reported by the technical staff: 
forwards (F), central midfielders (CM), central defenders 
(CD), wide players (WP): goalkeepers were excluded 
from the study. According to MEV performed in the ex- 
ercise test, players were divided into two groups (Hi-M: 
high-performers: able to run at a speed greater than me- 
dian value of all subjects; Lo-M: low-performers: able to 
run at a speed lower than that median). Furthermore, 
the subjects were divided into two groups according to 
median VEpeak (Hi- VE and Lo- VE, respectively). 

Linear regression analyses between VE and VO2, HR 
and I/O2, and HR and MEV, MEV and VO2 were com- 
puted at peak exercise. 

Two way analysis of variance (ANOVA) was carried out 
as follows: 1) dependent variable: VE peaky source of vari- 
ation: role, MEV; 2) dependent variable: FEVi; source of 
variation: role, MEV. All pairwise multiple comparison 
procedures were carried out using Holm-Sidak method. 

Students t-test was used for analysing statistical signifi- 
cance of differences in the following parameters: V02peak> 
Body Mass Index (BMI), peak Respiratory Rate {RRpeak)> 
peak tidal volume {Vtpeak)> MEV and V02peak/HRpeak on 
Hi-M vs Lo-M and Hi-VE vs Lo-VE, respectively. 

Factor analysis was computed on the set of cardiovas- 
cular, metabolic, and respiratory variables measured at 
peak exercise using non-rotated Principal Component 
Analysis (PCA). PCA is a simple, non-parametric method 
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of extracting relevant information from multivariate data- 
sets. The central idea of PCA is to reduce the dimensional- 
ity of a dataset consisting of a large number of interrelated 
variables, while retaining the variation present in the 
dataset. This is achieved by transforming into a new 
variable set, the Principal Components (PCs) which are 
uncorrected, and which are ordered so that the first 
few retain most of the variation present in all the original 
variables [37]. 

Multiple linear regression among variables resulting 
from PCA analysis was carried out. Correlation analysis 
between variables of interest was computed using Pearson 
coefficient. 

Statistical analysis was carried out using SigmaStat 
version 3.5 (Systat Software, Inc., USA), except for factor 
analysis carried out using SPSS vlO.l (SPSS Inc. USA). 

Results 

Anthropometric, demographic, and resting physiological 
characteristics of the whole study population and accord- 
ing to exercise capacity are shown in Table 1. The MEV 
median value was 18.65 km/h. BMI, height, weight, and 
age did not show any significant differences between 
Hi-M and Lo-M. MVY and FEVi were significantly higher 
in Hi-M than in Lo-M. There was no significant difference 
in anthropometric and demographic characteristics among 
different roles played. 

Table 2 shows the physiological parameters at peak 
exercise according to exercise capacity. Only VEpeai< was 
significantly higher in Hi-M than Lo-M, resulting from 
non significantly greater Vtpeai< and RRpeaio and non sig- 
nificantly lower Vd/Vt in Hi-M than Lo-M subjects. 

During the test, all players reached 97.9 ± 4.7% of their 
predicted maximal HR without any significant difference 
between groups. 

Table 3 shows resting and at peak exercise physiological 
characteristics according to VEpeak^ Median VEpeak was 
153.05 L/min. Hi-VE showed significantly higher 

V02peak, V02peak/HRpeak, MEV, BRR%, Vtpeak and FEVi 

than Lo-VE subjects. 



Table 1 Anthropometric, demographic, and resting 
physiological characteristics of the whole study 
population and according to MEV values 





Ail 


Hi-M n = 45 


Lo-M n 


= 45 


p 


BMI [Kg/m^] 


23.8 ±1.2 


23.6 ± 1 .4 


24.0 ± 


:0.9 


0.1326 


Height [cm] 


181.8 ±5.2 


184.0 ±4.9 


182.6H 


b5.2 


0.059 


Weight [Kg] 


79.4 ±5.5 


80.0 ±5.6 


78.8 ± 


:5.3 


0.3221 


Age [years] 


25.9 ±4.0 


25.2 ±3.5 


26.6 ± 


:4.4 


0.0969 


MW [L/m] 


201.1 ±21.7 


206.6 ± 20.4 


195.6± 


:21.6 


0.0155 


FEVi [L] 


5.0 ±0.5 


5.2 ± 0.5 


4.9 ± 


0.5 


0.0152 



Values expressed are shown as mean ± SD. Abbreviations as reported in the text. 



Table 2 Physiological parameters at peak exercise of 



subjects according to MEV 





All 


Hi-M n 


= 45 


Lo-M n = 45 


P 


MEV [Icm/h] 


18.5d 


il.l 


19.4± 


:0.6 


1 7.6 ± 0.7 


0.0000 


VEpeak[l-/min] 


151.7 d 


1 19.4 


156.21 


:19.8 


147.1 ±17.9 


0.0249 


(Vd/Vt)peak 


0.1 0± 


0.03 


0.09 ± 


0.04 


0.11 ±0.03 


0.0548 


RRpeak[breaths/min] 


53.3 d 


1 4.9 


54.51 


:5.1 


52.6 ±4.6 


0.0773 


Vtpeak[l-/min] 


2.87 ± 


0.45 


2.92 ± 


0.42 


2.84 ± 0.46 


0.3848 


HRpeak[b/min] 


186.0: 


±9.5 


186.2: 


t9.6 


185.7 ±9.4 


0.791 1 


V02peak[mL/min/Kg] 


63.3 d 


1 5.3 


63.31 


:5.3 


63.2 ±5.4 


0.9359 


BRR[L] 


41.9 ± 


20.8 


42.8 ± 


20.1 


41.0 ±21.4 


0.6726 


BRR% 


79.0 d 


1 9.3 


79.1 1 


:9.5 


78.8 ±9.1 


0.8750 



Values are shown as mean ± SD. Abbreviations as reported in the text. 



VEpeak was significantly correlated to V02peak (r = 0.619, 
p = 0.001) (Figure 1) and MEV was significantly correlated 
to V02peak (r = 0.267, p = 0.011). No statistical significant 
correlation was found between HRp^ak arid MEV. 

A weak, although significant correlation between 
Vd/Vt and VE at peak exercise was found in all players 
(r = -0.282; p = 0.007). No significant differences in V 
02peak > MEV and BRR were found among different 
roles, whereas mean MVV was significantly lower in F 
than in CM (193.3 ± 21.2 vs 206.1 ± 19.9 liters respectively, 
p = 0.043). ANOVA showed that the difference in the 
mean value of VEpeak and FEVi among the different levels 
of role of players was not significantly different after 
allowing for the effects of differences in MEV (p = 0.377 
and p = 0.543, respectively). 

Factor analysis using the PCA on the variables measured 
at peak exercise and at rest was carried out. Table 4 



Table 3 Resting and at peak exercise physiological 
characteristics of subjects according to l/Epea^values 





Hi-VE n = 45 


Lo-VE n = 45 


P 


MEV [km/h] 


18.9 ±0.9 


18.1 ±1.2 


0.0008 


VEpeak[L/min] 


1 66.3 ± 10.4 


137.1 ± 14.8 


0.0000 


(Vd/Vt)peak 


0.09 ± 0.03 


0.10 ±0.03 


0.0741 


RRpeak[breaths/min] 


53.8 ±5.1 


52.7 ±4.6 


0.2617 


Vtpeak[L/min] 


3.11 ±0.31 


2.62 ± 0.34 


0.0000 


HRpeak[b/min] 


185.8 ±8.9 


186.1 ±10.0 


0.8946 


V02peak/HRpeak[mL/min/bb] 


28.3 ± 4.0 


25.4 ±3.8 


0.0004 


V02peak[mL/min/Kg] 


65.7 ± 4.4 


60.8 ±5.3 


0.0000 


BRR[L] 


34.5 ±18.0 


49.3 ± 20.8 


0.0006 


BRR% 


83.4 ±7.8 


74.6 ± 8.7 


0.0000 


RRpeakW CO2/VO2] 


1.08 ±0.1 7 


1.10±0.10 


0.4576 


MW[L/min] 


208.0 ±18.0 


194.2 ±22.9 


0.0020 


FEVi[L] 


5.2 ±0.4 


4.8 ± 0.6 


0.0021 



Values are shown as mean ± SD. Abbreviations as reported in the text. 
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Figure 1 Relationship between minute ventilation and oxygen 
uptake at peak exercise as shown by the equation: VEpeak = 
41 .76 + (21 .88 * VO^peak) r = 0.61 9,p< 0.001 . 



reports the distribution along the first three components 
of the following variables: VEpeak > y02peak> (K//Vt)peai<, 

^Rpeah BRRpealo RRpealo V 02peak / ^Rpeak ^ FEVi. The first 

three components accounted for 74.33% of cumulative 
variance. Multiple linear regression on most representative 
variables on each component deriving from PCA at peak 
exercise is represented by the following equation: 

V02peak = 1.748 + (0.012 X VEpeak) 

+ (0.299 X FEVi)-{0m3 x BRR) R 
= 0.641 p < 0.05 (2) 

Equation (2) represents a statistically significant rela- 
tionship between V02peak and VE^eah EEVi and BRR. 

Discussion 

The novelty of this study is represented by the analysis of 
ventilatory parameters as evaluation tool of professional 
soccer players performance. Our elite soccer players were 
evaluated during seasonal activity. The main finding of 
our study is that VEpeak and FEVi but not V02peak are 

Table 4 PCA component matrix (peak values) 

Component 



1 



V02peak 


0.950 


0.014 


0.081 


VOipeak 


0.859 


-0.217 


-0.088 


VOipeak/HRpeak 


0.853 


-0.182 


0.305 


VEpeak 


0.720 


0.147 


-0.518 


RRpeak 


-0.305 


0.362 


-0.665 


HRpeak 


-0.180 


0.406 


-0.304 


BRR 


-0.197 


0.675 


0.688 


(Vd/Vt)peak 


-0.341 


-0.545 


0.301 


FEVi 


0.432 


0.817 


0.214 



the main determinants for discriminating high and low 
performers: nevertheless, neither Hi-M nor Lo-M were 
ventilatory limited. Therefore, our results are in agreement 
with available literature about elite athletes: the main 
limitation to maximal performance is represented by 
cardiovascular limit. Furthermore, our data highlight the 
key role of resting dynamic ventilatory parameters like 
FEVi as predictive factor of maximal exercise capacity, 
indicating the need to include lung function test in the 
routine evaluation of these elite soccer players. 

In our study population measured V02peak is in agree- 
ment with average values reported in previous studies on 
international elite soccer players [11]. Currently, ventilatory 
parameters are poorly investigated in these athletes and 
our study provides a contribution to this field. 

Significant relationships between VEpeak > FEVi and 
BRR with exercise capacity were found, and these pa- 
rameters were significantly different between high and 
low performers. This finding offers a new insight in this 
field as it provides a quantitative relationship between 
dynamic ventilatory parameters and exercise performance 
in these athletes. Indeed, V02peak described by a linear 
combination of most representative ventilatory parameters 
resulting by PCA analysis (i.e., VEpeak > FEVi and BRR) 
corresponds to a physiological relationship among these 
variables. Therefore, the physiological parameters obtained 
from a simple test like spirometry might be used to distin- 
guish high from low performers. 

The significant correlation found between {VJVt)peak 
and VEpeak on the whole study population represents an 
interesting finding: high values of VEpeak found in Hi-M 
elite soccer athletes correspond to low values of (K// V^)peai<« 
High ventilatory efficiency is more important in achieving 
better performance than VE. These parameters depend 
also on anthropometric and demographic characteristics 
which in our study were not significantly different between 
groups of players, independently of exercise capacity, 
ventilatory capacity, or role played. 

Our study did not find any significant difference in 
physiological parameters among different playing roles. 
This may be related to lack of specific training among 
the roles. 

Limitation of the study 

The subjective reason to stop exercise, either dyspnoea, 
muscular fatigue or both, was not recorded, despite most 
subjects informally reported muscular fatigue. This is 
not surprising given the lack of ventilatory limitation to 
exercise shown by these athletes. 

Conclusions 

Respiratory parameters can play a determinant role in 
qualitative and quantitative evaluation of professional 
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soccer players performance as it was confirmed by results 
discussed so far. Further studies should be addressed to 
investigate the impact of different training programmes on 
dynamic ventilatory parameters. 
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